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ABSTRACT 

We propose two new methods for measuring tangential peculiar velocities of rich clusters of galaxies. Our 
first method is based on weak gravitational lensing and takes advantage of the differing images of background 
galaxies caused by moving and stationary gravitational potentials. Our second method is based on measuring 
relative frequency shifts between multiple images of a single strongly lensed background galaxy. We illustrate 
this method using the example of galaxy cluster CL 0024+1654. 

Subject headings: galaxies: clusters: general- gravitational lensing 


1. INTRODUCTION 

The distribution of peculiar velocities of galaxies and clus¬ 
ters of galaxies is sensitive to the overall matter density in the 
Universe on large scales, and can be used to study structure 
formation. Peculiar velocities of galaxies have been used to 
reconstruct matter density helds (Colberg et al. 2000; Dekel 
& Ostriker 1999; Peacock 1999). Peculiar velocities of clusters 
could be used to determine large-scale bulk motion (Kashlinsky 
& Atrio-Barandela 2000) and check the reconstructed matter 
density helds. Direct methods to measure peculiar velocities of 
clusters of galaxies suggested so far are based on the frequency 
changes that moving clusters generate in the background pho¬ 
tons of the CMBR due to inverse Compton scattering or gravi¬ 
tational lensing. 

Measurements of cluster radial peculiar velocities are based 
on inverse Compton scattering of CMBR photons off hot elec¬ 
trons in the intra-cluster gas, the Sunyaev-Zel’dovich effect (Sun- 
yaev & Zel’dovich 1980). On average the CMBR photons gain 
energy from this scattering and a static and stationary cluster 
induces a decrement in the Rayleigh-Jeans part of the spectrum 
(the static thermal SZ effect). If the cluster is moving relative 
to the CMBR, an additional kinematic SZ effect is generated 
because of the bulk motion of the intra-cluster gas. The ratio of 
the kinematic to static thermal SZ effect, ATksz/ tSJssz, in the 
Rayleigh-Jeans spectral region may be expressed as 
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where is the radial peculiar velocity, L. is the electron temper¬ 
ature, and k is the Boltzmann constant, so that in moderate- or 
high-temperature clusters of galaxies the static effect will dom¬ 
inate. The kinematic SZ effect can be separated from the static 
effect using its different frequency distribution, even though 
their spatial distribution is the same (for recent reviews see 
Birkinshaw 1999 and Rephaeli 1995). The accuracy to which 
this radial velocity effect can be measured in a cosmological 
context was analyzed by Aghanim, Gorski & Puget (2001). At 
present only upper limits exist on radial velocities of clusters 
based on the kinematic SZ effect. Holzapfel et al. (1997) de¬ 
termined limits of peculiar radial velocities for two clusters, 
A2163 and A1689, with iltj ranges: -390 km s“' < Vj < 
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1860 km s“' and -460 km s“' < Vr < 985 km s“^. 

Tangential peculiar velocities of clusters can be determined 
using the moving cluster effect (Birkinshaw & Gull 1983; see 
also Pyne & Birkinshaw 1993 and Gurvits & Mitrofanov 1986). 
This effect is due to the fact that gravitational bending of light 
deforms the dipolar CMBR temperature anisotropy in the rest 
frame of the cluster, and thus when this pattern is transformed 
back to the CMBR rest frame (the observer’s frame) we do not 
recover the original, uniform, CMBR distribution. This effect 
is only of order 5 pK, far less than the static or kinematic SZ 
effects from prominent clusters (400 pK and 50 pK). Mea¬ 
surements of such small fluctuations in the CMBR on the few 
arcminute scale of clusters is beyond our current capability, and 
so no attempt has yet been made to measure this effect. Sensi¬ 
tivity improvements of a factor 100 are needed if this effect 
is to be measured, and this will be difficult in the presence of 
primary CMBR structures with the same spectrum. 

In this letter, we propose two new methods for measuring 
tangential peculiar velocities of rich clusters of galaxies with¬ 
out recourse to CMBR observations. Our first method is based 
on weak gravitational lensing. The distortions of images of 
background objects caused by a moving gravitational lens are 
different from those caused by a static lens. We calculate the re¬ 
sulting weak lensing signal and discuss how this effect could be 
used to measure tangential velocities of clusters using current 
technology. Our second method is based on frequency shifts 
generated by moving gravitational lenses, but instead of using 
the CMBR as a reference, we suggest using multiple images 
of a single strongly lensed background galaxy. We illustrate 
this second method using the example of galaxy cluster CL 
0024+1654. 


2. THEORY 

In the following derivation we use the thin lens approxima¬ 
tion for gravitational lensing and assume that the deflection an¬ 
gle is small (<C 1 radian). This approximation is well justi- 
hed for gravitational lensing caused by most clusters of galax¬ 
ies (see for example Schneider, Ehlers & Ealco 1992). In this 
approximation, the lens equation becomes 

6» = 6>^ + a, (2) 

where 9^ and 9 are the source and image coordinates, and a. is 
the observer’s frame bend angle (Schneider et al. 1992). We 
choose coordinate axes x and y to lie in the plane of the sky. In 
the rest frame of the cluster, the bending angle is 
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^ix,y) = {6^, 6y) = Vtp(x,y), (3) 

where x and y are angular coordinates in the sky and 'ip(x,y) is 
the 2-dimensional gravitational deflection potential. The change 
in the photon four-momentum in the observer’s frame is 


A = 




7'" = T>7" 


(4) 


where C, TZ, X and X>, are the four-dimensional Lorentz trans¬ 
formation, rotation (caused by gravitational lensing), identity 
and distortion matrices, 7 '" is the incoming photon four-momentum, 
and /3 = v/c is the peculiar velocity of the moving lens in units 
of the speed of light, c. 

First we briefly discuss the case when the gravitational lens 
is moving radially, i.e. I3 = vrIc. In this case, for small bending 
angles and peculiar velocities (5^1 and /3 <C 1), from Equa¬ 
tion (4), we obtain the deflection angle 


a={\+l3)5. ( 5 ) 

This result agrees with the results of Frittelli, Kling & Newman 
(2002; Equation 99). 

In theory this correction to the bending angle should be taken 
into account when using weak lensing to determine masses of 
clusters of galaxies. Cluster masses are determined from the 
bending angle. If the lens is moving radially, the bending angle 
changes because of its radial velocity, and the correct cluster 
mass, M, is related to the mass determined by ignoring the ra¬ 
dial motion of the lens, Mq, by 


M = 


Mo 

\+ (3 


(6) 


However, since we expect that the peculiar velocities of clusters 
are less than 3000 km s“', this effect is less than one percent, 
and is negligible relative to other errors even if we take it into 
account that linear structure formation theories underestimate 
peculiar velocities of clusters by about 40% (see Colberg et al. 
2000). Therefore we do not discuss this case further. 

We now turn to the effect caused by tangential peculiar ve¬ 
locities. For simplicity, we assume that the peculiar velocity of 
the gravitational lens, (3 = vt/ c is entirely in the plane of the 
sky. We choose x parallel to /3. The bending angle in the ob¬ 
server’s frame can be read from the spatial components of A 
(Equation 4, assuming <5^1 and (3 1) as 


“=(o ?)■*■ 

where 7=1 / a/I -/3^. This gives a. = 5 when (3 = 0, as ex¬ 
pected. The Jacobian, J, of the transformation from the source 
to the image plane can be seen from Equations (2) and (7) to be 


^ )-t7^2^(<^o), ( 10 ) 

72 -K + l\ J 

where the convergence, k, and shear, ( 71 , 72 ), are determined 
by the deflection potential, ip, and the second term, TZ^^((po), is 
a rotation of the source image by an unmeasurable orientation 
angle ipo, from the arbitrary assignment of axes x and y (see 
for example Peacock 1999). From Equations (9) and (10) we 
conclude that the shear is characterized by the two-dimensional 
vector held, 7 , with components 


ll = ^[-li^xx+tpyy] ( 11 ) 

J2=^[l1pxy + 'tpyx]- ( 12 ) 

Whenever these second derivatives of the deflection potential 
exist, we should be able to determine the tangential velocity 
(which enters here through 7 ) and the deflection potential sepa¬ 
rately since the potential is a scalar, continuous held. Unless the 
deflection potential is such that d^(V^ip)/dxdy vanishes iden¬ 
tically, 7 can be determined as 


l' 2 ,yy-Jl,xy 
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(13) 


where 71,7 and 72,7 are the double partial derivatives of 71 and 

72- 

The frequency shift in the spectrum of a lensed background 
galaxy caused by the gravitational held of a tangentially moving 
cluster of galaxies can be read from the time component of the 
four-vector A (Equation 4, again assuming (5 <C 1 and /7 ^ 1) 
as 


— =lf3Sx, (14) 

I'd 

where (3 = vt jc, and 5x is the x component of the rest frame 
gravitational bend angle (Equation 3). In the case of a gravita¬ 
tional potential with spherical symmetry, this reduces to 

Ai 2 / i/Q =-fPS cos (f, (15) 

where (p is the angle measured from the direction of the tan¬ 
gential velocity of the moving lens, which is identical to the 
expression derived by Birkinshaw & Gull (1983, their Equation 
9) and to first order in /J, to Equation (2) of Gurvits & Mitro¬ 
fanov (1986). 

Since there is a time delay between the two images of the 
strong lensed background galaxy, there is a small additional 
redshift difference, Az, between the lines of the two images due 
to the expansion of the Universe. This can be approximated as 


J = j2D_,_p2Z)^ (8) 

where is the two-dimensional identity matrix, and the two- 
dimensional matrix 


p2D ^ f Itpxx iPxy \ 
tpyx Ipyy J 


(9) 


(where 1/7 are second partial derivatives of the deflection po¬ 
tential with respect to variables i and j) describes the distortion. 
The distortion matrix can be decomposed as 


Az — Zlens) XI{ziens) At, ( 16 ) 

where H(ziens) is the Hubble parameter evaluated at the red- 
shift of the lens, ziens, and At is the time delay difference be¬ 
tween the two images. However, for a typical cluster, this red- 
shift, Az«5x 10 “" (l-l-z/cni) (//o/ 50 kms“' Mpc“')(Af/l yr), 
yields a velocity difference only of order 20 ms“' and is neg¬ 
ligibly small. The frequency shift in Equation ( 14 ) is not mea¬ 
surable unless one can find a reference frame, or measure dif¬ 
ferences in frequency shifts in multiple images of a single back¬ 
ground object. 
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Fig. 1.— Ax^ contours corresponding to 1, 2, and 3(T for the ellipti¬ 
cal potential of Equation (17). 2D cut: ellipticity, 77 and tangential peculiar 
velocity, V. The + sign shows the fiducial model parameters, 77 = 0.1 and 
V = 1500 km s~^ (see text for details of the simulation). 
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Fig. 2.— Ax^ contours corresponding to 1, 2, and 3it for an elliptical 
potential. 2D cut: normalization of the deflecting gravitational potential, ipQ 
and tangential peculiar velocity, V. The + sign shows the fiducial model pa¬ 
rameters, rpQ = 2.455 X 10“^ and V = 1500 km s“* (see text for details of the 
simulation). 


3. DETERMINING TANGENTIAL VELOCITIES 
3.1. Weak Lensing 

We carried out simulations to quantify the observability of 
the tangential peculiar motion on the shear field as revealed 
by images of background galaxies (Equations 11 and 12). The 
shear parameter, 7 , is directly related to the measurable ellip¬ 
ticity parameter, e = (ei,e 2 ), 7 ~ e, since k ^ 1. The ellipticity 
parameter as a function of position can be directly measured 
from the images of lensed galaxies using their isophotes (e.g., 
Schneider et al. 1992). We assumed that ellipticities could be 
measured for 40 galaxies per arcmin^, over an area of 50 x 50 
arcmin^, and that the ellipticity parameters of the background 
galaxies can be measured with 20% error. We average over 0.5 
arcmin^ cells, containing 20 galaxies each to obtain 5% ellip¬ 
ticity measurements in = 5000 samples. We choose an ellip¬ 
tical gravitational deflection potential, ip, which approximates 
the potential of a singular isothermal mass distribution at large 
radii (see for example Schneider and Bartelmann 1997) 


(17) 
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with 77 = 0.1, rcore = l^ and ipo = 2.455 x 10 “^ to show that the 
shape parameters of the potential (ipo, rcore, v) and the tangen¬ 
tial peculiar velocity of the cluster decouple, i.e. can be de¬ 
termined separately. We evaluated the errors using the Ax^ 
statistic, where we minimize 


i=N j=2 
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( 18 ) 


7 is calculated using Equations (11) and (13), and the elliptic¬ 
ity parameters, e are obtained by means of Monte Carlo sim¬ 
ulations assuming errors atj. S should be distributed like Ax^ 
with 2A^-4 degrees of freedom. Our simulations show that S is 
a non-degenerate four-dimensional ellipsoid near its minimum, 
so that P, Ipo, rcore, and p are independently determined, and 
that a tangential peculiar velocity of yj = 1500 kms“' can be 
measured to ± 300 kms“' (Icr, see Eigures 1 and 2). 


3.2. Frequency Shift 

The method of measuring the tangential peculiar velocity of 
a cluster of galaxies proposed by Birkinshaw Sc Gull (1983) and 
corrected by Gurvits Si Mitrofanov (1986) and Pyne & Birkin¬ 
shaw (1993) relied on observing the brightness change of the 
CMBR arising from the frequency effect. Equation (15). How¬ 
ever, it would also be possible to measure this frequency effect 
directly if there were a background with a sharp line feature. 
Such a background is available in the case of strongly-lensed 
multiple images of a single background galaxy. 

We illustrate this method using the well-studied, rich cluster 
of galaxies, CL 0024H-1654 (Shapiro & Iliev 2000; Tyson et 
al. 1998; Dressier, Gunn & Schneider 1985), at a redshift of 
0.395, which produces multiple images of a background galaxy 
located at a redshift of 1.675 (Broadhurst et al. 2000). The 
surface mass density of the cluster was found to be 
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with K = 7900 Mq pc“2, rcore = 35 kpc, and = 0.57 (Shapiro & 
Iliev 2000). The deflection angle for a spherical potential can 
be expressed as 


5(r) = 


AGM(r) 


( 20 ) 


where Ai is the mass included within cylindrical radius, r. The 
relative frequency shift between arcs in the direction of, and 
opposite to the direction of tangential motion of the cluster 
becomes twice as large as the single-image frequency shift in 
Equation (14). Assuming that the tangential velocity vt = 

1000 kms“E and is aligned with two of the lensed images, we 
obtain a frequency shift in emission (or absorption) lines equiv¬ 
alent to a Doppler shift of 1 kms“E To measure this effect we 
must find narrow lines from the same spatial location of the im¬ 
aged background galaxy in each arc, and measure the frequency 
difference between them. Narrow forbidden emission lines of 
[GUI] at AA4959, 5007, [Oil] at AA3727, 3729, and [Nil] at 
A6584, shifted into the infra-red (IR) band should be suitable. 

High signaEnoise spectra of such lines should be able to fix 
their central wavelengths to a few per cent of their widths. The 
precision of the measurement will be dependent on the stabil¬ 
ity of the spectrograph (which is better than 1 kms“' for stable 
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systems such as Phoenix on Gemini South; Hinkle et al. 2000) 
and the intrinsic widths of the lines being used. Narrow-line 
emitting regions within galaxies should have line widths of un¬ 
der 50 kms“\ and so spectra with signal/noise 10^ would be 
needed to fix the line centroids precisely enough to measure the 
cluster tangential velocity to 1000 kms“' or better. Such a high 
signal/noise cannot currently be achieved in the near IR at the 
magnitudes of interest for CL 0024-1-1654, but would be possi¬ 
ble with 50-m class telescopes. 

4. CONCLUSION 

In this letter we proposed two new methods for measuring 
tangential peculiar velocities of clusters of galaxies. Our first 
method is based on weak gravitational lensing. We carried out 
simulations to estimate how well we could determine tangen¬ 
tial velocities of clusters from distortions of the shapes of back¬ 
ground lensed galaxy images. We assumed that the ellipticities 
of background galaxies can be measured with 20% error. With 
today’s technology, obtaining such a large-scale and accurately 
measured shear field is possible only for a limited area. In order 
to measure this effect one would need to reduce the scatter in 
the intrinsic ellipticities of the background galaxies so that this 
noise source does not swamp the effect of tangential motion. 
One possible way of doing this would be to select the spherical 
components of lensed background galaxies, perhaps observing 
the bulge components in the IR for these galaxies using next 
generation space telescopes. 

Our second method is based on measuring frequency shifts 
between multiple images of a single strongly lensed background 
galaxy. This method allows us to measure components of tan¬ 
gential velocities of clusters in the direction of pairs of im¬ 
ages of strongly-lensed galaxies by measuring their relative red- 
shift. We used cluster CL 0024H-1654 to demonstrate that the 
needed velocity resolution is achievable with present-day spec¬ 
trographs. However, the signal/noise requirements are not achiev¬ 
able with existing near-IR instruments and telescopes, with im¬ 
provements of a factor 10 in throughput needed to make the 
measurement possible in this band. 

We might expect clusters of galaxies situated in superclus¬ 
ters to show the largest tangential velocities, since these clusters 
seem to move faster than "field clusters" (Colberg et al. 2000). 
While the differential redshift technique should still be feasible 
for such clusters, the shear technique may be compromised by 
the presence of the large-scale shear field of the supercluster. 

This work (SMM) was partly supported by NASA LTSA 
Grant NAG5-3432. 
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